Continuous mode cooling and phonon routers for 
phononic quantum networks 



S. J. M. HabrakenS K. Stannigel^ ^, M. D. Lukin^ P. Zoller^ ^ 
and P. RabF 

^Institute for Quantum Optics and Quantum Information of the Austrian Academy 

of Sciences, Innsbruck, Austria 

^Institute for Theoretical Physics, University of Innsbruck, Austria 
^Department of Physics, Harvard University, Cambridge MA, USA 

E-mail: Steven . habrakenQuibk . ac . at 



Abstract. We study the implementation of quantum state transfer protocols in 
phonon networks, where in analogy to optical networks, quantum information is 
transmitted through propagating phonons in extended mechanical resonator arrays 
or phonon waveguides. We describe how the problem of a non-vanishing thermal 
occupation of the phononic quantum channel can be overcome by implementing 
optomechanical multi- and continuous mode cooling schemes to create a 'cold' 
frequency window for transmitting quantum states. In addition, we discuss the 
implementation of phonon circulators and switchablc phonon routers, which rely on 
strong coherent optomechanical interactions only, and do not require strong magnetic 
fields or specific materials. Both techniques can be applied and adapted to various 
physical implementations, where phonons coupled to spin or charge based qubits are 
used for on-chip networking applications. 
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1. Introduction 

The successful application of laser cooling techniques for cooling isolated vibrational 
modes of micro- and nanomechanical devices [Il|2l[3llll[5l[6l[7l[8] has recently 
attracted a lot of interest in the control of macroscopic phononic degrees of freedom 
on a single quantum level. By now the preparation of mechanical resonators close to the 
quantum ground state has been achieved in different experimental settings [9l[T0l[TT] and 
coherent interfaces between mechanics and other quantum systems like superconducting 
qubits [9l [12] , spins [131 [E] or photons [15], |16] are currently developed. Beyond new 
possibilities to address fundamental questions in quantum physics [13 UHl [13 120], 
these experimental developments also provide the foundation for new, phonon-based 
quantum technologies. For example, in the context of quantum information processing 
and quantum communication, optomechanical (OM) slowing of light [211 ES] and first 
steps towards realizing a mechanical quantum memory [151 [IS| have been demonstrated, 
and the use of mechanical quantum transducers for interfacing different qubit systems 
[23] [211 125] has been proposed. For these applications mechanical systems benefit from 
the ability to interact with a wide range of other electric, magnetic and optical quantum 
systems, while still maintaining long coherence times and being compatible with scalable 
nano-fabrication techniques. 

The use of phonons for quantum information science applications is not new. 
Already in the first proposals for quantum computers, it has been suggested to employ 
vibrational modes of a trapped ion Wigner crystal for transmitting quantum information 
between spatially separated qubits [26j . More recently it has been shown that these ideas 
could equally well be applied in systems of coupled macroscopic mechanical resonators 
[23] [27] [28] , which extends the concept of a mechanical quantum bus to a wider range of 
atomic and solid-state systems. In analogy to optical fields, phonons can be confined in 
phonon cavities (e.g. represented by a high Q mechanical resonator), but also propagate 
freely along phononic waveguides. This suggests that many quantum communication 
and state-transfer protocols developed in the context of optical quantum networks 
[29l [30] l3T] could - on a smaller physical scale - also be implemented using acoustic 
phonons. Here, new approaches to design and pattern phonon waveguides based on 
phononic crystals structures [321 [33] provide a very promising and versatile platform 
for realizing such phonon networks in practice. However, compared with the relatively 
advanced field of optical quantum networks [30J, many equivalent control techniques 
still have to be developed for phononic quantum systems, which face the additional 
challenge that thermal noise in phononic channels is not negligible and would usually 
by far exceed quantum signals encoded in a single phononic excitation. 

In this work we address the problem of implementing quantum communication 
protocols in thermal phononic channels and show, how the addition of OM control 
elements can be used to realize a faithful transfer of quantum information between 
different nodes of the network. As a first key element to achieve this task, we describe 
the generalization of OM laser cooling techniques to multi- and continuous mode setups. 
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This approach is motivated by a recent proposal for interfacing OM systems (OMS) 
with phonon waveguides [25], and leads to a strong suppression of thermal noise within 
the relevant transmission bandwidth. Therefore, instead of pursuing the otherwise 
challenging task of cooling the whole network, this technique creates a 'cold' frequency 
window, which is sufficient to coherently transfer single quanta through an otherwise 
'hot' phononic channel. 

As a second control tool we describe the implementation of phonon circulators and 
switchable phonon routers, which enable a directed transfer of propagating phonons 
through large ID or 2D networks. In the optical domain circulators or other non- 
reciprocal devices are usually based on the Faraday-effect. In principle similar effects 
also exist for acoustic phonons in certain materials [M]. However, due to the required 
large magnetic fields and the use of materials with non-optimized mechanical properties, 
this approach is not suited for on-chip phonon quantum networks. Instead, we propose 
an integrated circulator for acoustic phonons that relies solely on an OM induced non- 
reciprocity [351 [36], where the directionality is imposed by the phase relation between two 
optical driving fields. Thereby, the device can be switched on or reversed conveniently 
and in combination with the above-mentioned cooling techniques, this coherent OM 
routing scheme is in principle sufficient to fully control the distribution of quantum 
information in large-scale phonon networks. 

The remainder of this paper is structured as follows. In section |2] we first present a 
brief introduction to phonon networks and the general input-output formalism, which 
is used to model these systems. In section |3] we revisit OM laser cooling and generalize 
it to multi- and continuous mode scenarios. As a basic application we then discuss in 
section |4] how an OM noise filter can be applied to transmit a quantum state through 
a thermal channel. In section [5] we describe the realization of phonon circulators and 
routers using coherent OM interactions. Finally, in section [6] we outline several potential 
systems for implementing phonon networks and then summarize the main results and 
conclusions of this work in section [T] 

2. Phonon quantum networks 

Figure [T^ shows a schematic representation of a generic phonon quantum network, where 
individual nodes i = 1, . . . , N are connected via a mechanical quantum bus. In analogy 
to optical quantum networks, each node contains an isolated solid state two-level system 
('qubit') with quantum information encoded in states |0) and The qubits interact 
with localized phonon modes ('phonon cavities'), which are in turn weakly coupled to 
a common set of propagating modes of a phonon waveguide or coupled mechanical 
resonator array ('quantum channel'). The full Hamiltonian of this system is given by 



N 




(1) 
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Figure 1. a) Schematic representation of a generic phonon quantum network. At each 
node qubits interact with a localized phonon mode, which in turn is weakly coupled to a 
phononic channel represented by an extended array of coupled mechanical resonators 
or a ID phonon waveguide, b) A phonon quantum network consisting of a set of 
localized vibrational modes, which are side-coupled to a common phonon waveguide 
and mediate the interaction between qubits and the propagating phonon modes (For 
a specific implementation of such a setup see reference |25j ) . c) The localized phonons 
decay into left- and right propagating modes in the waveguide with corresponding rates 
7i and 7^. The and ^^^^ denote the incident and outgoing waveguide fields at 
each node. See main text for more details. 



where H^^ode -^channel describe the dynamics of the individual nodes and the phononic 
channel, respectively, and i^int accounts for the coupling between them. In the following, 
we assume that the Hamiltonian for the individual nodes takes the form {h= 1) 

Kodeit) = oomblk + + Kit) {a\h + , (2) 

where the cj^ are Pauli operators and hi the bosonic operators for the local mechanical 
modes of frequency Um- The third term describes the coupling between the qubits 
and the local mechanical modes. For each node the qubit frequency splitting Ag(t) 
and the qubit-resonator coupling \i{t) can be tuned independently by changing the 
qubit frequency or modulating the coupling with local control fields (see section 
section [6|. For ~ ^\ this allows for a controlled mapping of the qubit state onto a 
phonon superposition, which we will use in our discussion of the state transfer below. 
The implementation of qubit-resonator interactions as given in equation ([2]) has been 
described for various charge- and spin-based qubits in the literature and in section [6] we 
briefly summarize some of the most relevant systems. 
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2.1. Phonon channels 

The mechanical quantum channel which is used to communicate between the nodes can 
in general be represented by a chain of A^^ch coupled mechanical resonators, 

-f^channel = zJ 2^ + + 2 z2 ~ Xg+i)"^ . (3) 

e=i e=i 
Here xe and pe are the position and momentum operators, m is the effective mass, 
Uq the bare oscillation frequency and the spring constant k accounts for the nearest- 
neighbor coupling. For small A^^ch the coupled resonators form a discrete set of collective 
eigenmodes. In this case quantum state transfer and quantum information processing 
protocols between two qubits can be implemented by addressing only a single collective 
mode, as has been discussed in the context of trapped ion systems [26] or coupled 
nanomechanical [23] and optomechanical [2H] resonators. 

In this work we are primarily interested in the opposite regime of extended arrays, 
A^ch ^1- In this limit -^channel cau be represented by a dense set of plane wave 
modes, -^channel = J2q^q^q^qy whcrc = Sq^q', Uq IS the phouou dispersion 

relation, and, for a lattice spacing a, the momentum label q is restricted to the first 
Brillouin zone q G (— This scenario is realized, for example, in large arrays of 
coupled nanomechanical beams [37] or in phononic band gap structures [32], where 
each resonator i corresponds to a vibrational mode of a unit cell, typically of size 
a ~ 1 yum. In the continuum limit we can identify a frequency range Ao; away from 
the band edges in which the coupled resonator array exhibits an approximately linear 
dispersion cu^ ~ cDq + c\q\, where c is the effective speed of sound and Uq a frequency 
offset. In this case it is convenient to introduce the normalized displacement field 
^{z) = ^r{z) + where 

•^^w = \/f E^^'^^' '^^(^) = JtT.^'""'^'^^ (4) 

describe right- and left-moving mechanical excitations of the phononic channel. Under 
the assumption of a linear dispersion these fields obey [^fi{z,t), $1^(2', t')] = and 

'^n/Liz, t), J/, t')] ^ e-^^«(*-*')5 (t-t'Tiz- z')/c) . (5) 

By assuming that the frequency of the local mechanical modes Um also lies within this 
frequency range, the coupling between the local resonators and the waveguide modes 
can be approximated by 

(6) 



ifint - Y I y dz J25iz - Zi) (bl^z) + h^^iz) 



where L is the length of the waveguide and Zi are the positions of the nodes along 
the waveguide (see figure [T|d). Below we identify 7 as the total decay rate of the local 
resonator modes into the waveguide. Note that equations ([s]) and ^ are valid for times 
|t — t'l > Au~^ and distances \z — z'\ > a. An explicit and more detailed derivation of 



these results is given in Appendix A for the case of a simple coupled resonator chain. 
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2.2. Input- output relations 

Under the validity of equations (|5| and ^ and in the hmit where cOm and the bandwidth 
Au are large compared to the other characteristic frequency scales, we can eliminate 
the waveguide modes and use an input- output formalism [38] to describe the effective 
dynamics of the coupled nodes. Using ^ and (|6| and making a standard Born-Markov 
approximation, we can derive a set of coupled quantum Langevin equations (QLEs) [38] . 
For each node we obtain 

where 7 = 7_r + 7l is the total decay rate of phonons into the waveguide. For side- 
coupled phonon cavities we would usually have •jr = 7^, but we below we describe 
scenarios where the emission either to the left or to the right is effectively switched off. 
In equation ([T]) we have defined incoming fields (see figure ^) 

&fen = lim$L(2:i + e), = lim $^(2;^ - e), (8) 

which specify the left and right moving waveguide modes before they interact with bi. 
Similarly, we introduce the scattered outgoing fields 

btnt = lim ^L{z^ - e), b^^^, = lim ^r{z, + e) . (9) 

The dynamics of the whole network is then described by the QLEs ^ together with 
the input-output-relations at each node, 

btntit) = bf.^it) + ^b,{t), (10) 

and the free propagation 

where Xj j+i = \zi+i — Zi\/c. Under realistic conditions we must also account for intrinsic 
phonon losses, backscattering and rethermalization, which we address below. 

2.3. Optomechanical control techniques for phonon networks 

The phonon network described above is formally equivalent to optical quantum networks 
and can in principle be used in a similar way to transfer quantum states from one node 
to another [2^ 1^ . However, in contrast to high frequency optical fields, the thermal 
population A^th = {Qy^v{.^m / ksT) — 1)^^ of the phonon modes is usually not negligible. 
For phonon frequencies in the MHz to GHz regime and temperatures T ~ 0.1 — 1 K 
the number of thermal noise phonons in the waveguide will exceed the quantum signal 
encoded in a single phonon. In addition, as shown in figure [l}:;, phonons emitted from a 
single side-coupled node will naturally propagate along two directions, i.e. 7^ = 7/j. In 
larger networks or 2D arrangements the inability to route propagating phonons severely 
limits an efficient implementation of quantum communication protocols. 

In the remainder of the paper we consider an extended setup as shown in figure |2] 
and describe how the integration of additional OM control elements can be used to 
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Figure 2. Basic setup illustrating the application of OM control techniques for 
implementing a directed quantum state transfer between two nodes through a thermal 
phonon network, (i) An OM noise filter is used to suppress incident thermal noise 
~ A'th within a relevant transmission bandwidth, (ii) OM phonon routers ensure 
a directed emission and complete reabsorption of individual phonons at the second 
node. To achieve the reverse state transfer the directionality of the routers can be 
reversed as described in more detail in section [5l 



overcome these fundamental limitations of phononic networks. The two key ingredients 
in this setup are: 

(i) An OM noise filter, which is used to suppress thermal noise within the relevant 
transmission bandwidth. 

(ii) Coherent OM phonon circulators or phonon routers, which allow for directed 
propagation and switchable routing of phonons through the network. 

3. Optomechanical noise filters 

In this section we describe the apphcation of OM laser cooling techniques to suppress 
thermal noise in an extended phononic quantum channel. Obviously, ground-state 
cooling of the whole network becomes inefficient as the system size increases, but it 
is also unnecessary, since usually only a few modes within a small bandwidth are used 
for transmitting quantum states. In the following we consider a scenario as shown in 
figure [3} where a single laser-cooled mechanical resonator provides a 'cold sink', and 
within a certain frequency range suppresses the thermal noise of the reflected modes of 
a continuous waveguide (or the transmitted modes in the case of a side-coupled phonon 
cavity). A similar configuration has been previously proposed for realizing a traveling- 
wave photon to phonon converter, where under 'impedance-matched' conditions an 
incoming photon is converted into a traveling phonon and vice versa |25]. In this sense 
cooling of the waveguide can be interpreted as a mapping of the optical vacuum onto 
the phonon channel, while in turn mechanical noise is upconverted into optical photons 
leaving the system. 
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Figure 3. Setup for an OM noise filter (OM continuous mode cooling). A mechanical 
resonator is coupled to a continuous phonon waveguide (7) and laser cooled via an 
optical cavity mode. Within a bandwidth ^ 7 around the mechanical frequency, the 
incident thermal noise can be efficiently suppressed, creating a dip in the filtered noise 
spectrum Np(uj) of the reflected waveguide modes. 



3.1. Single mode cooling 

Before addressing the suppression of thermal noise in extended phonon waveguides, let 
us first briefly review with the standard scenario for OM cooling [HUl EO], where the 
frequency Uc of an optical cavity mode is modulated by the displacement of a single 
mechanical mode. The optical cavity is driven by a coherent laser field of frequency 
Ud = oJc + ^ and in a frame rotating with the Hamiltonian for the OMS is given by 

-f^om = -5a^a + cumb^b + ga)a{b + b^) + iS{a) - a) . (12) 

Here the bosonic operators a and b represent the optical and the mechanical modes, 
respectively, ujm is the mechanical oscillation frequency, g the OM coupling constant and 
E the strength of the external driving field. Terms rotating with the optical frequency 
have been neglected by a rotating wave approximation (RWA). The OM interaction, as 



described by the third term in equation (12), derives from a linear dependence of the 



optical resonance frequency on the position quadrature of the mechanical mode [HI |32] ■ 
Including dissipation through cavity decay and intrinsic mechanical losses, the dynamics 
of the OMS is described by the QLEs 

d = {i5 — K)a — iga{b + 6^) + £^ — \/2Kain, (13) 

b = {-iujm - lo/'^)b - iga^a - ^To &o,in, (14) 

where k is the optical field decay rate and 70 = Wm/Qo is the mechanical damping 
rate for an intrinsic mechanical quality factor Qq. The 5-correlated operators and 
bo;m represent the vacuum input noise of the optical field and the thermal mechanical 
noise respectively. Within the relevant frequency range the latter is characterized by a 
non-vanishing equilibrium occupation number (&o,in(^)^o,in(^')) = NthS(t — t'). 

In the limit of strong driving the external field displaces both the optical and 
the mechanical modes by a large classical amplitude a = (a) = — i6{a)) and 
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(3 = {h) = —g\a^ juJrn-, where b{a) = 5 + 2g^\QS^ jujra- For |a| ^ 1, we can make a 
unitary transformation a — > a + a and b P + b and linearize the OM couphng around 
the classical mean values, 

Horn — —Sa^a + ujrnb^b + g{aa^ + a*a){b + b^) , (15) 

where we have redefined 6 (a) — )■ 6. If the external driving field is (near) resonant with 
the red mechanical sideband of the optical cavity, i.e. if 5 ~ — Wm, and if Um ^ i^, \gc(\ 
we can make a rotating wave approximation (RWA) with respect to and obtain a 
beam-splitter Hamiltonian 

Horn ^ i^bs = {uJm + Spb + g{aba^ + a*b^a) , (16) 

which describes a (near) resonant conversion of phonons into photons (and vice versa). 
Combined with the decay from the optical cavity, it allows for cooling: incident (noise) 
phonons are up-converted to photons and decay from the optical cavity. 

The QLEs for the linearized OMS can be conveniently solved in the Fourier domain 
as outlined in more detail in [Appendix B In the relevant weak coupling and sideband 
resolved regime {\g<y\ < k <^ Um = —S) we then obtain the standard result for the 
mechanical fiuctuation spectrum, 

(70 + lop)N 



{b\u)b{uj')) 



5{uj - uj') , (17) 



{U - UJmY + (70 + 7op) V4. 

where 7op = 2\ga\^ / k is the optically induced mechanical damping rate and N = 
-^th7o/(7o + 7op) + ^^/(4co'm) the reduced steady state resonator occupation number. In 
the following we are interested in OMS where both 7op ^ 70 and single-mode ground- 
state coohng N < 1 can be achieved. 



3.2. Optomechanical cooling in a multimode system 

The single mode cooling scheme can be generalized to cool mechanical networks 
consisting of a few coupled resonators only one of which is OM cooled, as shown in 
figure |4}d. In this case the linearized OM Hamiltonian is given by 

N 

Hom ^ -5a^a + ^ujrrfi]bj-^Kij {bib] + b^^b^^+ga (a)bi + afoj) , (18) 

where the Kij denote the nearest neighbor couplings and we have already performed a 
RWA assuming that Ki^j,g\a\, k <C oJm- The resulting QLEs can be written as 

d = {i5 — K)a — iagbi — v^Om (19) 

h = - {iuJm + lo/'2)bi+iKi^2b2-wa*a- ^/lobojri (20) 

bi<j<N = - {i(^m + 7o/2)6j + i{Kj^ijbj^i + Kjjj^ibj+i) - A/To&o'^in (21) 

bN = - {i(^m + 7o/2)&Ar + iKN-i,NbN^i - VTo&S , (22) 

where the fe^/^j^ are mutually uncorrelated noise operators associated with intrinsic 
mechanical damping of each mode. 



10 




Figure 4. a) Setup for OM cooling of a single mechanical mode and b) its 
generalization to multiple coupled resonators, c) The steady state fluctuation spectrum 
S{uj) as defined in the text is plotted for the 610 mode in an array of iV = 10 resonators, 
assuming that the 61 mode is OM cooled. The dashed line shows the result in the 
absence of OM cooling, i.e. g\a\ = 0, while the solid line corresponds to g\a\/K ~ 0.5. 
The other parameters for this plot are 6 = —uim, I0/K = 0.05 and k/K = 0.5. 



In figure |4p we plot the correlation spectrum (6jy(a;)67v(a;')) = S{uj)S{uj — u') for a 
chain of = 10 resonators with = K assuming that the bi mode is optically 

cooled. For reference the dashed line shows the undamped case ga = 0. In this 
case the heights of the peaks at frequencies Un = Um — 2K cos{mi / {N + 1)) are given 
by S{oOn) = (4iVth/7o)|cn(iV)|', where c„(j) = ^2/{N + l)sm{nj7T/{N + 1)) are the 
normalized mode distributions. When the OM coupling is turned on the peaks are 
broadened and their height is reduced, which corresponds to an overall cooling of the 
individual eigenmodes. However, we see that cooling occurs in a highly non-uniform way 
and only the few modes close to the cavity resonance are cooled efficiently. While the 
details depend on the ratios between K, g\a\ and k, we find that this behavior is quite 
generic and a precursor to the features we identify in the following for the continuous 
waveguide limit. 

3. 3. Optomechanical cooling of a mode continuum 

Starting from the multi-mode setting shown in figure |4]d, let us now address the limit 
— )■ cxD of a continuous mode waveguide as depicted in figure |3] by retaining the OM 
coupling to the first mode b = bi, but describing the other phonon modes 6j>i in terms 
of continuous left- and right-propagating fields ^lji{z). If K12 ^ -ft'i.i+i = K, we can 
adapt the input-output formalism introduced in section |2] to model the coupling between 
the localized mode b and the rest of the waveguide in terms of the scattering relation 

boutit) = 6in(t) + VlKt)- (23) 

Here &in(i) = — )■ 0,t) and &out(^) = ^r{z 0,t) are the incident and refiected 
waveguide fields and 7 ~ KI2/K is the characteristic phonon decay rate into the 
waveguide. Ignoring other, intrinsic loss channels for the moment, we obtain 

a = {i5 — K)a — iag{b + b^) — ^j2Ka\a (24) 
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Figure 5. Detailed results for OM noise filtering, a) Filtered phononic noise spectrum 
Np{lu) (black lines) and the corresponding photonic output spectrum (red lines), 



with the beam-splitter Hamiltonian ( f 6 ) (solid lines) and with the full linearized 
Hamiltonian ( 15 1 (dashed lines). The parameters are uJm/-f ~ 1200, 70 = 0, K/7 = 300, 
g\a\ = \/{"f + jo)k/2 and Nth — 40. b) Dip of the filtered noise spectrum Np{uj) for 
different values of the intrinsic mechanical damping rate 70/7 = 0, 0.0025, 0.005. 



b = — (icOm + 7/2)^ — ig{aa^ + a*a) — ^J^h^r^ , (25) 

for the linearized OM dynamics of the local mode and we see that the problem of cooling 
a continuous waveguide is formally identical to the single-mode cooling considered above. 
However, instead of looking at the occupation of the local mode b we are now interested 
in the spectrum of the reflected waveguide field 6out(^), given a thermal incident field 
(^L(^)^in(^')) = ^th^(^ — ^')- To do so, wc solvc the QLEs in Fourier space and write the 
result as 

Xut(u;) = S{u)Ai^{u) , (26) 
where Aj = (aj(w), aj(— w), 6j(u;), 6j(— cj))""" with j = "in", "out". The matrix S{u) is 



a 4 X 4 scat tering matrix and a more detailed derivation of equation (26) is given in 
Appendix B For given input noise correlations 6*^(0;, u') = (Ai^^u) ^ A]^{ijj')) we obtain 



the output correlation matrix Cout (1^5 f^') = (^out('^)®^out('^')) = S{u)Cia{u,u')S {lu'). 

3.3.1. Optomechanical noise filters. The quantity that we are interested in is the 
filtered noise spectrum Nf{uj) of the reflected field, which is defined by 

{blMboutico')) = Np{u)5{cu - uj') . (27) 

Figure |5^ shows the typical frequency dependence of Np{u). As expected, we see a 
strong suppression of thermal noise around the mechanical frequency Um- For 6 = —Um 
and under the condition \ga\ < k <^ Um, the dynamics of the OMS is well approximated 



by the beam-splitter Hamiltonian (16) and we obtain 



Nf(uj) ~ Nth { 1-^77 ^ 7 -t: ; -r 1 . (28) 



We see that, under these approximations, the 'impedance matching' condition 7op = 7 
[25] results in a complete cancellation Npiu = Um) = of the reflected noise on 
resonance. In this case the optical decay rate matches the mechanical waveguide 
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coupling and the OMS acts effectively as a double-sided phonon cavity with the thermal 
mechanical bath on one side and the optical T = bath on the other side. 



Under realistic conditions the noise cancellation described by equation (28) will 
not be perfect and to account for various imperfections in the system we assume in our 
discussion below a spectrum of the form 

Npiu) = iVth - (iVth - No)-, , (29) 

where 7 and cDm now include small corrections of the width and the center of the dip 
and Nq is a non-vanishing noise floor. Intrinsic limitations that lead to a finite A^^q 
are already contained in the full OM interaction itself, where for any finite Um off- 



resonant Stokes-scattering terms in the linearized Hamiltonian (15) induce corrections 
to the ideal beam-splitter interaction. However, as shown in figure [Sjo, by accounting 
for these effects up to second order in l/um we only obtain a small shift of the spectral 
dip Cjm — <^m — 5'^|«P/(2wm)^ and an offset A^o — l5'Qfp/(4c(;^), which is negligible 
under the weak coupling conditions mentioned above. A more crucial limitation for the 
OM noise filter arises from intrinsic mechanical losses 70 of the localized mechanical 
mode, which can be accounted for by introducing an additional decay rate 70 and 



associated noise operator 6o,in(i':) in the QLE (25) (see Appendix B). This leads to a 
finite A'q — 4A^th77o/ (7op + 7 + 7o)^- This means that, while for the optimized case with 
7 = 7o + 7op the local mode b can be highly excited, also the condition 7oA^th/7op <^ 1 
must be fulfilled to achieve good noise filtering. 

3.3.2. Propagation losses and scattering. The spectrum Np{u) determines the 
mechanical noise of the outgoing displacement field boutit) = ^r{z — )■ 0,t) immediately 
after the OM device and for an ideal phonon waveguide this noise spectrum will be 
the same at all positions z > 0. However, in a realistic setting, scattering losses in the 
waveguide and the associated noise lead to a rethermalization of the noise spectrum and 
Nf{uj) — 7- Nf{oo,z). For an approximately linear dispersion relation the propagation of 
the outgoing field can be modeled by 

^ + c|- ) ^r{z, t) = - {iQo + ) ^r{z, t) - ^<l'th(^, t) , (30) 



we outline a 



dt ' dzj ' ' V 2/^/ A 
where is the phonon mean free path in the waveguide and (^th{z^t) is a thermal noise 
field with {^\^{z,t)^tY.{z' ,t')) = Nt^6{z - z')5{t - t'). In [Appendix A.l 



microscopic derivation of this result for a coupled mechanical resonator array, where 
= c/70 can be connected to the intrinsic damping rate 70 of the individual resonators 
in the array. For 2; > and in a frame rotating with a)o we obtain 

^R{z,t)=e 2<.$^(^0,t--j -— = j dz'e , 



(31) 

and from this result we can derive the full noise spectrum along the waveguide 

Nf{u, z) = e-^Npiu, 0) + A^th (1 - e"^) . (32) 
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This result describes the rethermahzation of the noise dip on a length scale given by 
the phonon mean free path. For the relevant distances z <^ the noise spectrum 
can still be described by the Lorentzian shape given in equation (29), but the noise 
floor No{z) ~ Ai'o + {z/ly)Nth increases linearly with the distance to the OM noise 
filter. Note that similar quantitative conclusions follow from the standard approach 
of treating propagation losses in terms of a series of beam splitters [251 EH], but the 
present analysis also provides a direct connection to the underlying microscopic phonon 
scattering mechanism. 

To the extend that the scatterers are linear, phonon backscattering gives rise to 
small overall losses, which are on equal footing with the thermalization discussed above. 



4. Quantum state transfer in a thermal phonon network 

We now return to our original goal of implementing a quantum state transfer protocol 
between two qubits via an extended and thermally occupied phononic channel. For a 
simplified discussion we consider in the following the state transfer between two nodes 
assuming a unidirectional coupling, where 7 = •jr, 7l = and &i,in/out = ^fm/out 
shown in figure |2j In section |5] below, we describe how this condition can be realized 
using phonon routers. 



4.1. A tunable qubit network 

As a first step we describe the implementation of an effective qubit network with tunable 
qubit-waveguide couplings ri 2(t) by eliminating the dynamics of the local phonon 
modes. We start from the full set of QLEs, which, for each node j = 1, 2, derive 
from the Hamiltonian in ^ and are given by 

di = - iAg{t)ai + i\jit)aibj , (33) 
&i = -2zA,(t)K6,-6y:), (34) 

bj = - iumbj - -bj - i\j{t)a^_ - ^fybj^i^ . (35) 

The input-output relations are 

hj,ont{t) = &j,in(t) + Vlbj{t) and &2,in(i) = &i,out(i - 7-12) , (36) 

where 6i,in(^) = &in(^) describes the incident waveguide field before interacting with the 
first node. Hereafter, we absorb the retardation time T12 by redefining time-shifted 
operators and control fields for the second node j3H] and for simplicity we focus on the 
resonant case in which A^(t) = Um- 

We are interested in the regime in which the decay 7 of the the local phonon modes 
into the waveguide is fast compared to the coupling X{t). This allows us to adiabatically 
eliminate the phonon modes and to derive an effective description in terms of the qubit 



degrees of freedom only. In the frame rotating with ojm and to lowest order in A, (35) 
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can be solved to give 



T J — oo 



After reinserting this expression into equations (33) and (34) we obtain 

r,(t) 



yj_ ^ _ IIW 



r,(t)a^5,,i„ (3J 



&i = - r,(t) (1 + a^) + + 4i,al) . (39) 

Here Tj{t) = 4A|(t)/7 are tunable qubit decay rates and 

/t 
rf,e-i{t-)e-™«5^..^(5) ^ (40) 
-oo 

the associated effective noise operators which on the timescale obey 
[i?j^in(t), -Bjin(^')] ~ ^{t ~ t')- Provided that A(t) and cr^{t) also vary slowly on this 
time scale, we find that 

52,in(t) ^ -5i,in(t) + ^/^\(^)ai(^) . (4i) 

Thus we obtain a new set of effective quantum network equations for the qubit operators 
with tunable decay rates Tj{t). We emphasize that while in the following this effective 
description allows a simplified discussion of the state transfer protocol, it is not necessary 
to consider the regime A <^ 7 and a perfect state transfer can also be achieved using 
the full model ESI. 



4-2. Quantum state transfer through a phonon channel 

Our goal is to implement a quantum state transfer between the two qubits, i.e. 

li^o) = («|0)i + /3|l)i)|0)2 ^ |0)i(a|0)2 + /3|1)2) , (42) 

which is achieved via coherent emission and reabsorption of a single phonon in the 
waveguide. A solution to this problem has been first described for optical networks [29] , 
where it has been shown that by identifying appropriate pulses Tj(t), a perfect quantum 
state transfer can be implemented. Let us briefly summarize the main idea behind 
this scheme, for the moment assuming zero temperature. In this case the waveguide 
is initially in the vacuum state |vac) and the dynamics of the whole system can be 
restricted to the zero- and one-excitation subspace. Then, for an initial two qubit state 
l^po), we can deflne the amplitudes Vj(t) = (vac|(Oi, 02|cri(t)|'?/'o)|vac). From the QLEs 



(38) and (39) and the input-output relation (41), it follows that these amplitudes evolve 



according to 



Mt)= -^v,{t), (43) 

Mt)= -^vi it) - VUt)r2{t)v, (t) , (44) 
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and for initial amplitudes Vj{tQ), the general solutions is given by 

viit) = gi{t,to)viito) , (45) 
V2it) = Q2it,to)v2ito) + r{t,to)vi{to) . (46) 
Here Gjit^to) = e'^^^^^^^^/^ 

r(t,to) = - / dt' g2{t,t')^r^{t')r2{t')gi{t',to) , (47) 

J to 

is the state-transfer amplitude. For the initial state given in equation (42), ^2(^0) = 
and therefore a perfect transfer is achieved if at some final time tf the amplitudes 
approach 

gi{tf,to)^0 and \r{tf,to)\ ^ 1 . (48) 

While the first condition is easily fulfilled for sufficiently long, but otherwise arbitrary 
pulses, satisfying the second one requires control over their shape. Perfect state transfer 
is only possible if the total number of excitations is conserved, i.e., if no population gets 
lost from the one-excitation subspace. This means that the norm gf{t, tQ) + T'^{t, to) = 1 
for all times, which after taking the time-derivative of this condition yields 

Vf;(t)v,{t) + v/fVt)i;2(t) = . (49) 

This result also follows from the requirement that the total scattered field after the 
second node vanishes at all times, i.e. -B2,out(t)|V'o)|vac) = 0. Therefore equation (49) is 
often referred to as the dark-state condition. 

A set of pulses Ti{t) and r2(t) that realize a perfect state transfer can always 
be found numerically by choosing Ti(t) such that ^i(t/,to) and then solve 
equations (43) and (44) iteratively, choosing r2(t) at each time such that the dark- 
state condition (49) is fulfilled. Further, simple analytical expressions for pulses that 



realize a perfect state transfer may be obtained from a time-inversion argument 
Without loss of generality, we can assume that [tf,to] = [—Tp/2,Tp/2], where Tp 
is the pulse length. A control pulse Ti{t) for the first qubit gives rise to a wave 
packet a{t) = ^yTl{t)Ql{t, to)vi{to) in the waveguide. For reasons of time-inversion 
symmetry, the inverted wave packet a{—t) is fully absorbed if the reverse pulse ri(— t) 
is applied. This implies that, in the special case of a time-symmetric wave packet 
a{t) = a{—t), the wave packet generated at the first node will be fully absorbed if 
we choose T2{t) = ri{—t). Describing the symmetric wave packet by the differential 
equation a{t) = g{t)a{t), where g(t) is anti-symmetric git) = —g{—t), one can derive 
the following differential equation for the pulse-shape Ti{t) from (43) [15] : 

Fi = + 2g{t)T^ . (50) 

For the simplest choice of g(t), i.e. g(t) = sign(t)(rmax/2), where Fmax is the maximal 
decay rate, the solution reads 



r,{t) = { \2-e-rn...tJ^m.. li t < T^it) = r,{-t) . (51) 

Fmax ff t>0 

We will use this specific pulse shape for our numerical examples discussed below. 
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Figure 6. Quantum state transfer through a noisy channel, a) Spectral overlap of 
the transfer pulse |F(a;)p and the filtered noise spectrum Np{uj) for A'q = and 
Tmax/T = 0.10. b) Effective occupation number N^g as a function of 7/rniax for 
No/Nth — 0.05. c) Numerical simulation of a quantum state transfer in the presence 
of thermal noise with iVth = 0.5. d) The same state transfer with the OM noise filter 
switched on and Fmax/T = 0.10. In both plots the blue and red solid lines show the 
excited state populations pl{t) — {a^_^{t)a'^_{t)) for the qubits 1 and 2, respectively. 
The dotted lines show the corresponding pulse shapes Ti{t) — Ti[t) /T-ai_^y_ as defined in 



equation (51), but with cut-offs. In d) the dashed line shows the results obtained from 
a two qubit cascaded master equation with an effective thermal occupation number 



4-3. State transfer through a noisy quantum channel 

Now we consider the same state-transfer problem but in the case in which the in- 
field of the quantum channel is characterized by a non-vanishing noise spectrum 
{bI^{uj)Bi^{uj')) = N{u)6{u - u'), which can either be a flat thermal background 



A^(ci;) = Nth or the flltered noise dip N{uj) ^ Nf{uj) as discussed in section 3.3 



4-3.1. Effective thermal occupation number. In the presence of thermal excitations in 
the quantum channel we must necessarily go beyond the single excitation subspace and 
a closed analytic solution to the state transfer problem is no longer available. To gain 
more intuition on the impact of noise on the state-transfer fldelity, let us flrst consider 
a single qubit. We assume that at time t = to the qubit is prepared in its groundstate 
and we then switch on the coupling to the waveguide Ti{t), for example, using the 



pulse shape defined in equation (51). In the regime where the noise amplitude is low. 



N(uj) <^ 1, we can linearize the QLEs (38) and (39) and obtain for the final qubit 
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excitation — Ncs ^ 1, where 

iVeff= / dt,dt2 ^T^{tr)Tr{t2)Gi{tfM)Qi{tM){Bl{t^)B,^{t2)). (52) 

In the case of incident 5-correlated thermal noise and assuming that the pulse duration 
is sufficiently long, i.e. tf — to ^ r~!^^, we find Nes = A^th- This means that during 
the state transfer an erroneous excitation or de-excitation probability for the qubits 
of ~ A^th can be expectec|f] In turn, for the filtered noise spectrum Nf{uj) defined in 
equation (29) we obtain 



oo 



iV., = / du: \F{u:)\^NA^) ^ ^^^o + T N,^ ^ ^^^^ 

J —CO '^1 ~r ^ max 

where in general F(uj) = J^^ dt e"^"^^ ^/Ti(t)Gi{t f , t) and the integral in equation (53) 
has been evaluated for the pulse shape given in equation (51). This shows that using 
the OM noise filter to clean the waveguide, the effective occupation number can be 
considerable reduced compared to A^th, assuming that A^o 1 and that the bandwidth 
of the transfer pulse fits within the noise dip. This is illustrated in figure |6j where we 
plot the spectral overlap between |F(w)p and Npiuj) and the resulting A'efr for different 
ratios 7/rmax- 



4-3.2. Discussion. To verify our approximate analytic predictions we convert the 
cascaded QLEs into an equivalent cascaded master equation [38] and simulate the full 
state transfer numerically. Since a master equation for the two qubits can only be derived 
for 5-correlated noise, we include the OM noise filter as a third subsystem to emulate 
the spectral variations of the incident noise, as outlined in [Appendix C[ In figure [6|3 and 
figure [6]i we simulate the transfer of a single qubit excitation from node 1 to node 2 and 
compare the case of a thermal quantum channel where N{u)) = A^th to the case where 
the OM noise filter is switched on and N{uj) = Nf{uj)- We see that even for A'th = 0.5 
thermal noise in a phonon quantum network already almost completely washes out the 
transferred signal, while in combination with the noise filter a high fidelity state transfer 
is still possible. In figure [6[i we also compare the full numerical results with a master 
equation for the qubits only, assuming a (5-correlated noise with an effective occupation 
number A^eff- We find that apart from small corrections when pulses are switched on 
and off, this reduced model describes the state transfer very accurately, which shows 
that, indeed, A'^eff is the relevant parameter for a quantum channel with a filtered noise 
spectrum. 

As a specific example, we consider a potential realization of phonon networks using 
the on-chip OM structures discussed in reference [25]. We assume local phonon modes 
of frequency Um = 27r x 4 GHz and a decay into the waveguide of 7 = 27r x 25 MHz. 
For a quality factor Qo = 10^ the intrinsic decoherence rate is 70 = 27r x 4 kHz. For 
these parameters we plot in figure [7] the fidelity J-" for transferring the superposition 
state \ip) = (|0) + \l))/\/2 between two nodes of a phonon network as a function of 

I Of course, for larger A^th the non-linearity of the qubit must be taken into account. 
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Figure 7. Quantum state transfer fidelity T for a superposition state = 
(|0) + |l))/\/2 as a function of the thermal occupation A'th of the waveguide and 
different values of Fmax- For this plot we have used the pulse shape given in equation 



(51 1 and assumed 70/7 = 1-6 x 10 



A^th and for different effective qubit decay rates 70 <^ Tmax < 7- For the numerical 
simulations we have used the effective model with A^^efr as obtained in equation (53) 
and with Ni^ = •yo^th/l- The fidelity is defined as the overlap between the target state 
and the actual state after the transfer, i.e. J-" = Tr{ptarP(^/)}- We see that compared 
to the case where no cooling is applied and the fidelity already drops significantly for 
A^th ~ 0.1, the presence of the OM noise filter can push this bound to much larger 
occupation numbers. In particular, in the example considered here, this means that 
instead of requiring temperatures of T < 100 mK, the OM cooling scheme allows the 
implementation of high fidelity state transfer protocols at much more convenient liquid 
Helium temperatures T = 4 K where A^th ~ 20. Note that for the assumed value of 
Tmax = 10~^7 = 271 X 250 kHz, the total transfer time of ~ 1/is is compatible with 
decoherence times that are achievable with various different solid-state qubits. A few 
specific examples will be discussed below in section [6] 



5. Phonon routers 



In the previous sections we have studied the transfer of single excitations through a 
thermal phonon network assuming that two nodes are coupled in a unidirectional way. 
This situation is applicable only to specific setups, for example, when the two nodes are 
located at the two ends of a single waveguide. In general the emission of phonons into 
left and right propagating modes, reflection and interference effects, or imperfect routing 
of phonons through multi-port junctions in a 2D setting will limit the implementation 
of efficient state transfer protocols in larger networks. In optical networks, many of 
these problems can be overcome by using optical circulators and optical isolators for the 
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Figure 8. A phonon circulator, a) General scheme for realizing a three port circulator 
based on three coupled phonon cavities, where one of the tunneling amplitudes is 
complex. A cyclic scattering of phonons between the three ports is achieved for t — 
and (/3 = 7r/2. b) A two mode OM setup for implementing an effective non-reciprocal 
tunneling amplitude te^'f . The two coupled optical cavities mediate an effective hopping 
between the phonon modes bi and 62, where the overall phase is contorolled by the 
phases 0i,2 of the external driving fields, c) Scattering probabilities \Si^j[(jj)\^ into 
the different ports j = 1, 2, 3 of the circulator for an incoming signal of frequency lo in 
port 1. The solid line shows the results for an ideal device (70 = 0) and the dashed 
lines the results for a finite intrinsic loss rate 70 = 7/20 for each of the three local 
phonon modes. 

directional routing of photons. In the context of OMS it has aheady been suggested 
to use the intrinsic non-hnearity of OM interactions to induce non-reciprocal effects for 
light [321 EH]- In the following we describe a related scheme, which allows us to engineer 
coherent non-reciprocal effects for phonons, where the directionality is simply imposed 
by the phase difference between two optical driving fields. 

5.1. A three-port phonon circulator 

To illustrate the basic idea of a non-reciprocal phonon device let us consider the 
minimal setting of a three port circulator as shown in figure |8^. The localized phonon 
modes 6^=1,2,3 are mutually tunnel-coupled with a strength \t\ and coupled to phonon 
waveguides with a decay rate 7. The Hamiltonian is given by 

i^circ = ^-^1^* + ^ (^i^^e'^ + hh\ + hh\ + H.c.) , (54) 
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assuming that the resonator frequencies Um are all equal. The crucial term in this setup 
is a complex tunneling amplitude te^'^ between two of the resonators, which cannot be 
absorbed into local redefinitions of the bi. Therefore, this phase can be associated with 
an effective magnetic field for the phonon modes. Previously, such a setting has been 
described for superconducting microwave cavities, where the effective magnetic field 
can be engineered using superconducting qubits and external fluxes [lU |15] . Below we 
describe a non-magnetic approach to achieve this complex hopping amplitude in OMS. 

Including the coupling to the waveguides and in a frame rotating with the resonator 
frequencies Um, the QLEs read 
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(55) 



with the input-output relations &j,out = i>j,m + ^/lb,j. The above set of QLEs can be solved 
in Fourier space and the input-relation can be applied to obtain the full scattering matrix 
S{uj) of this system. By choosing (p = n/2 and t = 7/2, such that the decays into the 
waveguides are matched to the tunneling amplitudes, we find that for frequencies around 
u ~ Um (w ~ in the rotating frame), the scattering matrix is given by 



(56) 
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Up to factors i, which can absorbed in the definitions of the operators, Eq. (56 ) describes 
a perfect circulator. For (p = — 7r/2 we obtain a similar result, but with the scattering 
direction reversed, i.e. 6i,out = ^&3,in, &2,out = &i,in and b^^out = ib2,m- In figure ISb we 



plot the scattering probabilities \Si^j{u)\ , defined as &j,out(i^) = <Si^j{uj)bi^in{uj), as a 
function of the frequency u (relative to Um) and for t = 7/2. We see the emergence 
of the ideal circulator relations close to resonance. We also find that these features are 
robust when we add an additional intrinsic loss rate 70 = 7/2O. 



5.2. Optomechanically engineered non-reciprocity for phonons 

To implement a complex tunnel amplitude te'^'^ between two localized mechanical modes 
we propose a setup as shown in figure [8]d. Here, the localized phonon mode 63 is coupled 
to modes 61 and 62 mechanically with a (real) tunneling amplitude t, while the tunneling 
between bi and 62 is mediated by two driven optical cavities. The Hamiltonian for this 
system is 

3 

Hcirc = J2 ^'"^^^^ + ^^^1^3 + ^2&3 + H.C.) + + H^ra ■ (57) 
1=1 

In the frame rotating with the frequency of an external driving field the Hamiltonian 
for the two coupled optical cavities is 

ifc = ^ -^idldi - J(ai4 + aia2) + ^ ^ {ci\e"''' - H.c. j , (58) 

i=l,2 i=l,2 
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where 8i are the strengths and the phases of the optical driving fields. Finally, 

^^om = 5Z + (59) 

i=l,2 

describes the local OM interactions. 

As discussed in section |3} in the limit of strong driving we can change to a displaced 
representation and linearize the OM coupling around the classical expectation values 
Oil = {tti). In the present case the classical field amplitudes are given by 

^ (k - i52)£ie"f'' + iJ£2e'*' 

{K-i5i)£2e"t'^ +iJ£ie"''' 
"2 = —f ^TT • . N I j2 ^ (61) 

where we have assumed that both cavities decay with the same rate rate k and we have 
absorbed non-linear shifts of the detunings by a redefinition of 5j. We write = \ai\e'^'^' 
and obtain the linearized OM coupling 

Hom = Yl 9\c^^\ie-''''a^ + e'^^a}^{h + &!)• (62) 

i=l,2 

Note that for given 5j, J and n the external control parameters £i and 0i provide enough 
fiexibility to independently adjust the phases ipi and keep |ai| ~ \a2\. 

We focus on the specific case 5i = 5, \ai\ = a and introduce symmetric and anti- 
symmetric optical modes a± = (ai ± a2)/\/2 with detunings 5± = 5 ± J. Further, we 
assume that —S± ~ Um g\c(i\, which allows us to make a rotating wave approximation 
with respect to Um- Then, after changing to a frame rotating with Um, the total 
Hamiltonian is given by 

^circ = -J2 + ^(^1^3 + ^2&3 + H.C.) + 

where Aj^ = 6^ + u)m- In a final step we assume \A,y\ ^ ga and treat the OM 
coupling using second-order perturbation theory. Apart from small frequency shifts 
for the mechanical modes this results in an effective tunneling Hamiltonian 

Hcirc ^ tbibl + th2h\ + teffe*'^6i4 + H.c. , (64) 

where = ipi — ip2 and 

Thus, by choosing the external control parameters such that = ±7r/2 and that tes 
matches the bare mechanical tunneling amplitude t, the setup shown in figure [8)d realizes 
a switchable three-port phonon circulator as discussed in the previous subsection. 
Note that the interaction with the optical modes also introduces an additional decay 
channel with rate 7op ~ g'^a'^n/A'^, where A"^ = A^^ + Al^. Compared to the noise 



[aV(e*'^i6i + e^^2 62) + al(e*'^^fei-e*'^2 62) + H.c. , (63) 
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filtering scheme described above, this requires slightly lower cavity decays satisfying 

K <^ \A\ <t: Um- 

As a specific example let us consider the phonon waveguide scenario discussed above 
with typical phonon frequencies Um ~ 27r x 4 GHz and a phonon-waveguide coupling of 
7 = 27r X 25 MHz. If we choose Si = —ujm we obtain t^s = 9^o? I J and the conditions 
^eff = 7/2 can be achieved for (70; ~ 27r x 110 MHz and J = 27r x 1 GHz. For the 
same parameters a value of k < 27r x 50 MHz is sufficient to suppress the optically 
induced decay rate 7op = 2g'^a'^K/.P below the value 7op/7 < 0.05 shown in figure [s]:. 
This is only slightly below demonstrated values for k in state of the art OM crystal 
cavities [2^125]. 

6. Implementations 

The fabrication and control of mechanical systems, resonator arrays and phonon 
waveguides as well as their coupling to other microscopic quantum systems (qubits) is 
currently a very active field of research. For many of these systems the general control 
techniques described in this work could provide the basis for phonon-based quantum 
communication applications or mechanical quantum interfaces in hybrid qubit setups. 
In the following we present a brief discussion of several potential candidate systems for 
implementing phonon networks. 

6.1. Phonon channels 

As described in section [2| a ID phonon channel can be realized by fabricating a large 
array of coupled nanomechanical resonators with a bare frequency Um and nearest 
neighbor coupling g < Um- This scenario has been experimentally studied, for example 
in reference where the resonators were charged up and coupled via electrostatic 
interactions. Alternatively, the resonators can also be coupled mechanically via bridges 
or the support [M]. Both approaches are suitable for realizing phonon channels with 
frequencies Um ranging from MHz to a few 100 MHz, where also very high Q- values 
around Q ~ 10^ — 10^ are nowadays routinely achieved. At T ^ 100 mK this frequency 
range corresponds to a few tens to a few hundred thermal phonons, which can still be 
efficiently suppressed using the proposed OM noise filtering scheme. 

A more practical and very versatile approach for implementing phonon waveguides 
based on phononic bandgap materials has recently attracted increasing attention [321 133] . 
Here, a 2D structure with periodically varying mechanical properties is designed such 
that a complete band gap in the mechanical dispersion relation appears. Then, phonon 
wave guides can be realized by introducing line-defects in these structures, to which 
the phonons are confined as long as their frequency lies within the bandgap of the bulk 
phononic crystal. With this approach phonon waveguides with frequencies Um ~ 1 — 10 
GHz can be realized, where even at T = 4 K the thermal occupation of the waveguide 
is only a few tens of quanta. Further, as indicated in Figure [1]d, such waveguides can 
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be easily combined with localized phonon cavities and integrated OM devices [25] 



6.2. Qubit-phonon interfaces 

The implementation of phonon quantum networks also requires the realization of 
coherent and controllable interfaces between mechanics and stationary qubits. Here 
the ability of mechanical systems to respond to weak optical, electrical and magnetic 
forces enables the coupling of mechanical resonators to a large variety of spin or charge 
based qubits and a few prototype examples are outlined in the following. 



6.2.1. Spin qubits. Qubits encoded in localized spin degrees of freedom, for example 
in spins associated with Nitrogen vacancy (NV) impurities in diamond [17] or phosphor 
donors in silicon [l8], can be coupled to mechanical motion using magnetized tips 
[T^t HU HH]. A strong magnetic field gradient V-B from the tip induces a position 
dependent Zeeman shift of the spin resulting in an interaction of the form 

^nodc = ycr^ + (^rnb^b + A(6 + 6"^)cr^ + ^ cos(a;™^t)cr^.. (66) 

Here uq ~ GHz is the bare spin splitting and A = gsHB^o^B/ (2h) is the Zeeman shift 
per zero point motion oq, where Qg — 2 and /xb is the Bohr magneton. For nano-scale 
mechanical resonators with frequencies in the few MHz regime this coupling can be 
substantial and reach values A/ (27r) !^ 10 — 100 kHz [30] ■ Base on this coupling the use 
of mechanical transducers to mediate spin-spin interactions in small resonator arrays has 
been proposed previously [23], and the present techniques can be used to extend these 
ideas to larger networks. To achieve an effective Jaynes-Cummings type interaction as 
given in Eq. (|2]), the spin is driven by a near resonant microwave field with frequency 
= ojq + 5 and Rabi frequency VL as described by the last term in Eq. (66). Then, 
by changing into a dressed spin basis and making a rotating wave approximation with 
respect to ujq, the effective interaction is given by [ISj [50] 

i^node ^ + ojJ)% + Asin(^)(a+6 + a_6^), (67) 



where = ~ ojm and tan(6') = Vl/5. The qubit-resonator coupling can be 

controlled by adiabatically varying the effective detuning or the mixing angle 9. 

Instead of using external magnetic field gradients, alternative schemes to strongly 
couple spins to mechanical motion have been recently suggested for carbon nanotubes 
[^ [52] . Here a single electron is localized on the nanotube and couples to its 
vibration via spin-orbit interactions. In this case even larger couplings A/(27r) 0.5 
MHz and larger mechanical frequencies > 100 MHz are achievable. Although a 
controlled fabrication and positioning of carbon nanotubes is still challenging, a phonon 
quantum bus could be realized by electrically coupling the nanotube to other mechanical 
resonators, which can be fabricated in a more controlled manner. 
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6.2.2. Superconducting qubits. The coupling of nanomechanical systems to various 
types of superconducting qubits has been demonstrated in recent experiments [9llT2|l53]. 
While for superconducting qubits other ways for communication exist, for example 
via microwave transmission lines, the coupling to phononic channels could still be 
interesting for creating interfaces to other quantum systems, especially to optical qubits 
for long-distance quantum communication [2H [25]. Depending on the type of qubit 
('charge', 'phase', 'flux', ...) the qubit-resonator interaction can be due to electrostatic 
[T71 IM] . piezo-electric |9] or magnetic interactions |55l |56l |57], and can in all three 
cases be substantially larger than in the case of spin qubits. Since the bare transition 
frequency of superconducting qubits is typically in the GHz regime, a resonant coupling 
to mechanical systems in the MHz range can again be achieved as described above, 
namely by applying additional driving fields to compensate for the frequency mismatch 
[5H 158] . Alternatively, superconducting qubits could be coupled directly to high 
frequency mechanical modes as demonstrated in Ref. [9]. 



6.2.3. Quantum dots and defect centers. The bending of a mechanical resonator locally 
deforms the lattice configuration of the host material and by that shifts the electronic 
states associated with artificial quantum dots or naturally occurring defect centers in 
solids. In the bulk this deformation potential interaction is usually responsible for 
line broadening of optical transitions and other decoherence effects, but in the case of 
confined modes it can also lead to a significant coherent coupling to single phonons. In 
Ref. [59] the coupling of a quantum dot to the fundamental bending mode of a doubly 
clamped beam has been considered, leading to a deformation potential coupling of the 
form 

i/def = X{b + 6t)|e)(e|, hX ^ (H^ - (68) 

where |e) denotes the electronically excited state, / the beam length, t its thickness and 
Oo the zero point motion. The Sg and are deformation potential constants for the 
ground and excited electronic states and for typical values Sg ^ ~ 1 — 10 eV and micron 
sized beams a coupling strength of A/27r ~ 10 MHz can be achieved. This is already 
comparable to the radiative lifetime Fg of the electronically excited state, but can in 
principle be pushed to values A ~ 1 GHz using much smaller, phononic Bragg-cavities 
as suggested in Ref. [60] . 

As a potential scenario where the deformation potential interactions could be used 
to implement a controllable qubit phonon interface, we consider an NV center embedded 
in a diamond nanoresonator. The electronic ground state of the NV defect is a spin 
triplet, where two states |0) = \ms = —1) and |1) = \ms = 1) can be used to 
encode quantum information. The qubit states in the electronic ground state, which 
are highly immune against external perturbations, can be coupled via an optical Raman 
process involving an electronically excited state, which, in contrast, is strongly coupled 
to phonons [^I] . In Appendix D we adiabatically eliminate the dynamics of the excited 
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state and derive an effective spin-phonon interaction of the form 

i/eff ~Aeff(t)K& + fT_6t). (69) 

Here the a± are Pauh operators for the qubit subspace and for optimized laser detunings 
Ks(t) = 4:\Qo(t)Qi(t) / co^ is a tunable coupling where ^o,i{t) are the optical Rabi 
frequencies. Note that this coherent interaction is accompanied by dissipation at an 
average rate feff(t) = 4:Xflo{t)fli(t) / u'^ and the ratio Xcs/^cs = A/F is not affected by 
the off-resonant Raman coupling. This is in contrast to cavity QED where the optical 
field couples to the atomic coherence and not to the population of the excited state as 
described by i/def- 

7. Conclusions and outlook 

In summary we have described the implementation of dissipative as well as coherent OM 
control elements for realizing quantum communication protocols in extended phononic 
networks. We have shown how OM continuous mode cooling schemes can be used to 
create a cold frequency window in the thermal noise spectrum of a phononic channel and 
we have analyzed the fidelity of quantum state transfer protocols under these conditions. 
Further, we have proposed the realization of non-reciprocal phononic elements, which 
rely on strong coherent OM interactions and where the directionality is simply controlled 
by the phase of an external laser field. Based on this principle, various switches and 
routers for propagating phonons can be constructed, which allow for the implementation 
of efficient quantum communication protocols also in larger phonon networks. 

Both the OM noise filter as well as the phonon router can be realized with 
state-of-the-art mechanical systems. Combined with the ongoing developments in 
the control of qubit-resonator interactions, they could soon provide the elemental 
tools for implementing phononic quantum networks. As potential applications of 
such networks we envision the distribution of entanglement within a larger quantum 
computing architecture, where propagating phonons can replace or complement direct 
qubit shuttling techniques [621 163] . Here the ability of mechanical systems to interact 
with various different types of qubits, makes phononic channels particularly suited for 
hybrid qubit settings. Beyond quantum communication applications, the OM control 
techniques described here may also be used to probe the propagation and scattering of 
single-phonon wave packets. 
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Appendix A. Coupled resonator cirrays 

In this appendix we derive the effective propagation equations and input-output relations 
for phononic quantum channels consisting of a large array of A^ch ^ 1 coupled mechanical 
resonators. Assuming a homogenous system the resonator array is described by the 
Hamiltonian 

-f^channel = ( 2^ + 2^^^°^^ ) ^ 2 ^ ~ (^-1) 

£=1 ^ ^ e=i 

and the equations of motion for the position operators xe are given by 

xe = -iolxe - ujoK{2xe - xe+i - xe^i), (A.2) 

where K — k/{muJo) is the nearest neighbor phonon-phonon coupling strength. For a 
large array we can assume periodic boundary conditions and write 



^^(^) = 7^ E V ^ (e^^-^/^-6„(t) + e-^-^/^-6t (^)) , (A.3) 

where bn{t) — fc^e""^"* are bosonic operators for the plane wave modes labeled by 
n — — (Ach/2 — 1), . . . , A^ch/2 and normalized to [6„, foj^,] = Sn,n'- The eigenfrequencies 
are given by 



= Y'^o + "^Kuo [1 - cos(27rn/A^ch)] • (A.4) 

We are interested in the regime where the total length L as well as the other relevant 
scales of the network are large compared to the spacing a between the individual 
resonators. We introduce a continuous field 



such that X£ = = ai). The quasi- momentum q is restricted to the first Brillouin 



zone q E {—n/a,n/a] and coq = a/c^q + '^Kcuq [1 — cos(ga)]. 

To proceed we now consider the tight binding limit K -C coq, where Uq ~ 
Wo + K{1 — cos(ga)) and the phonon modes form a band of width Aou — 2K around a 
large center frequency ui — ujq -\- K . Further, for frequencies around ui the dispersion 
relation is approximately linear and can be written as 

(jjq'^^u + c{\q\ - vr/(2a)) = loq + c|g|, (A. 6) 

where c = Ka is the sound velocity and Cjq = ujq — (7r/2 — 1)K a frequency offset. 
Therefore, as long as we are interested in the dynamics of phonon modes away from the 
band edges we can set i?channei — X^gl'^o + c\q\)h'^qhq and approximate the displacement 
field by 

m^^[m + ^\z)]. (A.7) 
V Au; 
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Here xq = ^/h/2muj and the field ^{z) = a/2c/L e^'^^bq is normalized to 

(A. 



-iwoit-t') 



5 { At ] + 6 i At + — 



[^z,t),<i>\z',t')] ^ ^ ^ 

where At = t — t' and Az = z — z'. The field operators can be decomposed into a right- 
and a left-moving component $(2) = ^r{z) + ^Liz), as defined in equations ^ and 
@. 

The coupling of the localized phonon modes to the phonon channel can be written 

as 

1 ^ 

Hint = 2 ^^'^(^^ ~ ^^)^' ('^•^) 

where ki^e = k\oc is non-zero only for site £, which are next to a local mode i (we 
can assume side-coupled resonators, such there is only one neighbor). We define 
K\oc = fcioca^o^^o where xq the zero point motion of the local mode. Then, under the 
assumption that Um ~ w 3> -ft'ioc, we can make a RWA and obtain the coupling given 
in equation ^ with a resonator decay rate given by 7 = 2Kf^^/ Au. 



Appendix A.l. Propagation losses 

To model propagation losses in the phonon waveguide we add an intrinsic loss channel 
with rate 70 for each of the waveguide modes. Making the RWA right from the beginning, 
the dynamics of the whole resonator array is modeled by the coupled QLEs, 

b, = - (i{ujo + ir) + ^) 6, + ^ + - v^6^,in(t), (A.IO) 

where bg is the bosonic operator for the mechanical resonator at site i and the 
b£,in(t) are uncorrected thermal noise operators in(^)^f',in(^')) = NthSe/'6(t — t'). 
In the corresponding momentum representation bq = 1/ y/N^h e*''"^6^ and 6q_in = 
^*'"^^^,in) the coupling is diagonal 

= - (^Wg + y) &g - VTo&g.inW, (A.ll) 

where as above coq = uo + K{1 — cos(ga)) ^ ujq + c\q\. From this equation we obtain the 
propagation equation for the right moving field ^ji{z,t) = a/2c/L ^^^^ e'"^^bq{t)^ 

{j^+^Yz) "^^^-^'^^ = " ('"^o + y) '^«(^'^) - v^'^'th(^,t). (A.12) 
Here the thermal noise field is defined as 



^th{z,t) = 




(A.13) 



such that it is normalized {^l^{z,t)^th{z' ,t')) = Nti,6{z - z')5{t - t'). 
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Appendix B. Optomechanical cooling 



After linearizing tlie OM coupling the QLEs (13) and (14) can be written as 



dtA{t) = -MA{t) - V2/€An(t) - v^5o,in(t), (B.l) 

where we have grouped operators as A{t) = (a(t), a'l'(t), 6(t), Ai^(t) = 

{ain{t),4nit),0,0V and Bo,in{t) = {0,0,bo,in{t),bl -^{t)f . The matrix M is given by 

/ —i5 + K iga iga \ 

i6 + K, —iga* —iga* 



M 



\ 





iga* 
—iga 



iga 
-iga* 



iu +2^1 




We define the Fourier-transformed operators as a{uj) 
and obtain 



(B.2) 

[a(a;)]t 

A{uj) = -X{u) [V2kAUuj) + ^oBo,M) , (B.3) 

where we set X{uj) = [M. — iujl]~^ and now A{u) = {a{u),a^{—uj),b{cL!),b'^{—uj))'^, etc. 
The non-vanishing correlations of the noise operators are {ai^{ijj)al^{u')) = 6{u — u'), 
{bo,i^{uj)bl .^{uj')) = (Ath + l)6{uj - u') and (6j .Ju;)6o,in(c^')) = Nti,6{u - u'), within the 
relevant frequency range. Then, the stationary fluctuation spectrum of the mechanical 
mode is given by 

io-uj').{BA) 

Under the weak coupling and sideband resolved condition this expression simplifies to 



the result given in equation (17) 



In section 3.3 we are interested in spectrum of the scattered field &out(^) in the 
case where the OMS is in addition coupled to the phonon waveguide with rate 7. 
We define A{u) and the intrinsic noise Bo^ni^) as above, but we set A^=ia^out{^) = 
{au{u),al{—u),bu{u),bl{—ijj))'^. The input-output relations can then be written as 
^out('^) = Aia{u) + \/TZA{u) where TZ = diag(2K, 2k, 7, 7) is a diagonal matrix. Then, 
together with equation ( |B.3 ) we obtain 

A,M = Siuj)AUuj) - S'iuj)Bo,U^). (B.5) 

where S{ijj) = 1 — 7lX{u) and S'{ijj) = ^/yo}ZX{uJ) and in the definition of X{uj) we 
replaced 70 — > 7o + 7- The filtered noise spectrum of the out-field is given by 

Nf{lo) = |5(u;)32p+Ath(|<S(u;)33|'+|5(a;)34|'+|5'(a;)33|'+|5'(a;)34n, (B.6) 



where in equation (28) the approximation Nf{ui) ^ Ath|iS(u;)33| has been made. 



Appendix B.l. Multi-mode cooling 

We can use the same approach to solve for the stationary state of the multimode 



system described in section 3.2 In the following we assume for simplicity the validity 
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of the RWA, such that the QLEs for the equations for the annihilation operators 
A = (a, 61, . . . bjy)'^ form a closed set 



dtAit) = -MA{t) - V2KAUt) - v^So,in(t), (B.7) 
where now iin(t) = (ai„(t), 0, . . . , 0)^ and Bo,Ut) = (0, . . . , ^^'^^(t))^ and the 



matrix Ai can be derived from equations (19)-([22j). Following the same steps as above 
we obtain 

{blituMco')) = 7oiVth \Xitu),+,,+if^ 5{lo - tu'), (B.8) 
which we have used to evaluate the multi-mode cooling results presented in figure 
Appendix C. Cascaded master equation 

To simulate the state transfer between two nodes in the presence of incident noise and 



OM cooling, we map the qubit QLEs (38) and (39) with time dependent decay rates 



Fi(t) and F2(t) onto an equivalent cascaded master equation [3Sj, which we can then 
integrate numerically. Since we cannot treat the filtered input noise directly using a 
master equation we consider a unidirectional network where the OM cooled phonon 
cavity is included as a first system. The OM cooling is modeled by an additional decay 
channel for the phonon cavity at a rate = 2g'^\a\'^ / k, which is matched to the decay 
rate 7 of the phonon cavity into the waveguide. The cascaded master equation, which 
describes this system is given by 

p = -i[H, p] + (A^th + l)V[S]p + N,^V[S^]p + ioMh]p , (C.l) 

where T>[c\p = cpc^ — {c^cp + pc^c)/2. By identifying cq = &, Cj=i,2 = o"- and Fq = 7 
the collective operator S = Y.k=o,i,2 v^c^ and H = Y.k>i V^^Wi{clci - cjck). Ath 
is the thermal occupation number of the incident white noise. 



Appendix D. NV-phonon coupling 

In this appendix we show how for an NV center in diamond the deformation potential 



coupling can be used to realize the effective spin-phonon interaction (69). The NV 
center has S = 1 triplet ground state and we assume that the qubit is encoded in the 
states |0) = \ms = —1) and |1) = \ms = +1). The qubit states are coupled by external 
laser fields to an electronically excited state |e), which is coupled to the deformation of 
the local lattice structure induced by the vibration of the beam. In the frame rotating 
with the laser frequencies the Hamiltonian is given by 

H=J2 [A,|j)(j|+^-(t)(|j)(e| + |e)(j|)]+a;„6t5 + A|e)(e|(6 + 6t),(D.l) 
i=o,i 

where the ^j{t) are tunable Rabi frequencies and we have chosen the zero of energy to 
coincide with the excited state level such that Aq and Ai are detunings of the drive 
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lasers from the |0) to |e) and |1) to |e) transitions, respectively. The interaction term 
~ A can be eliminated by a polaron transformation 

H^e^He-^ with S = —\e){e\{h^ -h) , (D.2) 

which transforms the Hamiltonian into 

H = ujj)% + ^ A,|j)(j| + [(fio|0)(e| + ^^l|l)(e|) A + H.c] (D.3) 

i=o,i 

and is still exact. Since the ratio rj = X/um ^ 1, we can expand to first order in rj and 
obtain A ~ 1 — 77(6''' — b). 

Our goal is to induce a coherent Raman transition from state |0) to state |1) while 
simultaneously absorbing a phonon from the nanomechanical oscillator. To make this 
process resonant we set Aq = A + c<;.m/2 and Ai = A — Um/'2, where A = (Aq + Ai)/2 
is the overall detuning. We change into an interaction picture with respect to Hq = 
Umb'^b + X]j=o 1 write the total state as \4')(t) = X]j=o 1 e '^^^ 
equations of motion are then given by 

dtmit) = - y |^e)(t) - t (1 + ^^""'^^ - r/e-^-'"*6) \i>,){t) (D.4) 

i=o,i 

^t\^i'J=o,l){t) = - t^je'"'^' (1 - r]e'''-V + r^e-'^^'b) (D.5) 

where we have added an imaginary part to model the radiative decay of the exited 
state with rate Fg [§j Assuming that |Ao,i| ^ ^o,i) we can approximately integrate the 
equation of motion for \%lje){t), 

and insert this result back into the equations of motion for |'?/'o,i)(^)- By keeping only 
non-rotating terms, we obtain for the iV'o) subspace. 



d^it) = --^\i^o)it) 
Ao + 



(D.7) 



and a similar result for the evolution of |?/'i)(t). In the limit |Aj|, |Aj ± w^l ^ Tg and 
using the definition Ao,i = A iWm/Z, we can identify an effective spin-phonon coupling 
Aeff and effective decay rates F^^ for each spin level, 

A2-a;y4' " ^e^- (0-8) 

The ratio between the coherent coupling and the mean decay rate rcfr = ivf^ + r|,^^)/2 
is optimized for A = where Xes/^es = A/Fg. 

Here we have assumed that the NV center is coupled to a single vibrational mode. 
In general, nanostructured resonators will support multiple mechanical modes so that 

§ For simplicity we here omit the associated recychng terms, which can be included, for example, using 
a stochastic wavcfunction formalism [38] . 
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the phononic part in the Hamiltonian (D.l) generahzes to -ffphonon = Sfc^^fc^I^fe + 



\^){^\^k ^k{bk + bl). However, since the higher-frequency modes of nanostructured 
resonators are separated from the fundamental mode by ~ GHz, the other modes are 
highly off-resonant and their contributions to the resulting spin-phonon coupling are 
negligibly small. See reference ^2] for a similar discussion. 
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